Systematic study of three-nucleon force effects 
in the cross section of the deuteron-proton breakup at 130 MeV 
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High precision cross-section data of the deuteron-proton breakup reaction at 130 MeV are pre- 
sented for 72 kinematically complete configurations. The data cover a large region of the available 
phase space, divided into a systematic grid of kinematical variables. They are compared with the- 
oretical predictions, in which the full dynamics of the three-nucleon (37V) system is obtained in 
three different ways: realistic nucleon-nucleon (NN) potentials are combined with model 37V forces 
(3NF's) or with an effective 3NF resulting from explicit treatment of the A-isobar excitation. Alter- 
natively, the chiral perturbation theory approach is used at the next-to-next-to-leading order with all 
relevant NN and 37V contributions taken into account. The generated dynamics is then applied to 
calculate cross-section values by rigorous solution of the 37V Faddeev equations. The comparison of 
the calculated cross sections with the experimental data shows a clear preference for the predictions 
in which the 3NF's are included. The majority of the experimental data points is well reproduced 
by the theoretical predictions. The remaining discrepancies are investigated by inspecting cross 
sections integrated over certain kinematical variables. The procedure of global comparisons leads 
to establishing regularities in disagreements between the experimental data and the theoretically 
predicted values of the cross sections. They indicate deficiencies still present in the assumed models 
of the 37V system dynamics. 

PACS numbers: 21.45.+V, 25.10.+S, 21.30.-x, 13.75. Cs 



I. INTRODUCTION 

The dynamics of the three-nucleon (37V) system can be 
very accurately studied by means of the nucleon-deuteron 
breakup reaction. Its final state, constrained by only gen- 
eral conservation laws, provides a rich source of informa- 
tion to test the nuclear Hamiltonian. It is of particular 
importance for components of the models which account 
for subtle effects, like three-nucleon force (3NF) contri- 
butions to the potential energy of the 37V system. Precise 
predictions for observables in the 37V system can be ob- 
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tained via exact solutions of the 37V Faddeev equations 
for any nucleon-nucleon (TV TV) interaction, even with the 
inclusion of a 3NF model pj. To investigate details of 
the dynamics of the 37V system, in addition to elastic 
Nd scattering data, reliable deuteron breakup data sets, 
covering large regions of the available phase space, are 
needed. Unfortunately, it still remains difficult to per- 
form such measurements at the required level of preci- 
sion. In our previous paper we have started to report 
results of a project dedicated exactly towards such an 
aim. Here we continue with the presentation of a sys- 
tematic set of breakup cross-section values and we com- 
pare them with theoretical predictions based on various 
dynamical assumptions. 

Properties of fcw-nuclcon systems at not-too-high en- 
ergies are determined mainly by pairwise nucleon-nucleon 
interactions. Models of 7V7V forces describe the long 
range interaction part according to the meson-exchange 
picture, while the short range is based on phenomenol- 
ogy, and adjusted by fitting a certain number of parame- 
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ters to the NN scattering data. The present generation 
of NN potentials reaches an unprecedented accuracy in 
describing the pp and np observables below 350 MeV, 
expressed by a \ 2 P er degree of freedom very close to 
1. In fcw-nucleon studies the most widely used so-called 
"realistic" NN potentials are Argonne Vis (AV18) Q, 
charge dependent (CD) Bonn jjJ-Q- Nijmegen I and II 
(Nijm I, Nijm II) [6]. Their full equivalence with phase 
shift analysis guarantees that all two-body aspects of 
the interaction are taken into account when these NN 
force models are used in microscopic calculations of few- 
and many-nucleon systems. 

At the more fundamental level of Quantum Chromo- 
dynamics (QCD), the strong force between the nucle- 
ons is understood as residual color force. A direct de- 
scription of fcw-nucleon systems at low energy from first 
principles would require the solution of QCD in the non- 
perturbative regime which is not possible at present (ex- 
cept on the lattice). On the other hand, the low-energy 
dynamics of QCD can be studied in the chiral effec- 
tive field theory (EFT) framework. This is a system- 
atic approach which incorporates the spontaneously bro- 
ken approximate chiral symmetry of QCD and is based 
on the most general effective Lagrangian for Goldstone 
bosons (pions in the two-flavor sector of the up and 
down quarks) and matter fields (nucleons, A-resonances, 
...). In the pion and single-baryon sectors, S'-matrix 
elements can be calculated in chiral perturbation the- 
ory (ChPT) via an expansion in terms of {Q/A. x ) v - in 
powers v of a low-momentum scale Q, associated with 
small generic external momenta and with the pion (light- 
quark) mass. Here, small means with respect to the scale 
A x , corresponding to the chiral symmetry breaking scale 
of the order of 1 GeV. Motivated by successful appli- 
cations of ChPT in the irir and irN sectors, Weinberg 
proposed to extend the formalism to systems with two 
and more nucleons, where non-perturbative calculations 
are necessary to deal with the shallow bound states (or 
large scattering lengths) H,|i|. According to Weinberg, 
ChPT can be applied in that case not to the amplitude 
but to a kernel of the corresponding dynamical equa- 
tion which may be viewed as an effective nuclear poten- 
tial. Few-nucleon S'-matrix elements are generated non- 
perturbatively by iterating the potential in the dynami- 
cal equation. The first application of this approach in the 
2N sector was performed in [iH. ITlj. At present, the 2N 
system has been studied up to next-to-next-to-next-to- 
leading order (N 3 LO) in the chiral expansion 0, 
while the three- and more numerous nucleon systems 
have so far been analyzed up to next-to-next-to-leading 
order (NNLO) HI E^. It should be stressed that this 
approach offers an unique possibility to estimate uncer- 
tainties of the theoretically predicted physical quantities. 

High-quality models of the NN potentials, when ap- 
plied to calculate observables in the 3N system, revealed 
discrepancies between the pure pairwise dynamics and 
the experimental results. The most promising and widely 
investigated explanation is the presence of three-nucleon 



interactions. The realistic potentials are therefore sup- 
plemented by 3NF models, usually refined versions of the 
Fujita-Miyazawa force {vfy, m which one of the nucleons 
is excited into an intermediate A via a 2-7r-exchange be- 
tween both remaining nucleons. The most popular ver- 
sion of such an interaction is the Urbana IX [T^l force. 
The Tucson-Melbourne (TM) [l<| 3NF extends this pic- 
ture by allowing for additional processes contributing to 
the pion rescattering at the intermediate nucleon. An al- 
ternative mechanism of generating a 3NF is based on the 
so-called explicit A- isobar excitation poLl2lll22ll23| . Cal- 
culations are performed in a coupled-channel approach 
and the effective 3NF is generated (together with other 
A-isobar effects) due to the explicit treatment of the de- 
grees of freedom of a single A. Finally, within the ChPT 
framework both 2N and 3N forces (as well as nuclear 
currents) are derived from the same effective chiral La- 
grangian and are thus fully consistent with each other. 
This leads to a consistent model of the NN and 37V in- 
teractions, which also strongly constrains the parameters 
of the 3NF. As stated earlier, presently the results in 
the 3iV system are only available at NNLO. The analysis 
at N 3 LO requires sophisticated analytical and numerical 
calculations. This work is in progress. 

The role of 3NF effects has been recognized already 
when studying the bound states of three nucleons. No 
realistic potential approach can reproduce the binding 
energies of 3 He and 3 H when the calculations are based 
on NN forces only [24|. When 3NF contributions are 
taken into account, the 3 H and 3 He binding energies 
can be described accurately (by construction, because 
parameters of the 3NF are usually fitted to match the 
triton binding) - see e.g. [25j. These combined models 
of NN and 3N forces also describe the 4 He binding en- 
ergy, indicating that 4NF are presumably small j2(J . For 
the description of the level schemes of p-shell nuclei, the 
most simple 3NF's show failures which motivated more 
sophisticated 3NF models leading to encouraging agree- 
ment between theory and experiment |27|. Here, we will 
restrict ourselves only to the models mentioned above. In 
the isospin T = 1/2 state, they are expected to be very 
similar to the extended versions used in [27]]. An analo- 
gous conclusion is obtained within the ChPT framework 
- inclusion of 3NF graphs leads to an improved descrip- 
tion of few-nucleon bound states |l4|. Further evidence 
of relevant consequences originating from introducing ad- 
ditional dynamics into the 37V system comes from the 
coupled-channel approach - the binding energies of 3 He 
and 3 H are much closer to the experimental values when 
the A-isobar contributions are included and the differ- 
ence of the two bindings is well matched |2^| . 

Presently, the richest evidence for the importance of 
3NF effects is deduced from the elastic nucleon-deuteron 
scattering observables. The picture emerging from the 
comparisons of various data with theory is, however, 
rather ambiguous. In several cases where the TV TV forces 
alone fail to reproduce the observables the inclusion of 
3NF's leads to significant improvements psl 12^. I30L l3ll 
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El HI HI (for earlier references c.f. Q). Alas, 
in several cases discrepancies between the experimen- 
tal data and theoretical predictions remain, even if the 
presently available full 3TV dynamics is taken into ac- 
count. This statement is especially true for various po- 
larization observables pffll l3ll l33l l33 | , but holds also for 
certain cross-section angular distributions (see e.g. |34j). 
Those failures, confirmed by different calculational ap- 
proaches, indicate that the 3NF models are still missing 
some relevant ingredients, while for the ChPT framework 
they might suggest the necessity for including higher or- 
der (at least N 3 LO) terms for the 3TV system. 

Since the theoretical models clearly need more con- 
straints from the experimental data, it is natural to ex- 
tend the investigations of the 3TV system to the nucleon- 
deuteron (Nd) breakup reaction. The continuum of the 
final states, which has to be simultaneously described 
in its full richness by the assumed dynamical model of 
TV TV and 37V interactions, should provide a lot of infor- 
mation to pin down the details of the theoretical models. 
Unfortunately, this field has hardly been explored ex- 
perimentally and only at lower energies, below 30 MeV 
nucleon energy (see [lj and [3 71 fo r references; the most 
recent results can be found in |38j and |39|. and in |4p| | 
at much higher incident energy). In the region of inter- 
mediate energies (30 MeV - 100 MeV) only at 65 MeV 
several isolated kinematical configurations have been in- 
vestigated with respect to cross sections and analyz- 
ing powers 0, l42l l43l ISf . Comparison of those data 
with the theoretical pred ictions obtained within the ap- 
proaches m m |2Zl discussed above shows again a 
mixed picture: sometimes the agreement is improved by 
including 3NF's, in some cases the 3NF effect is negligible 
and there are cases in which inclusion of 3NF's moves the 
prediction away from the data. Since the thorough the- 
oretical study [45( of the full phase space of the breakup 
reaction shows that significant effects can be expected, 
there is a strong need for data which precisely and sys- 
tematically scan large ranges of the final state kinemati- 
cal variables. 

Therefore, we have performed a 1 H(d,pp)n breakup ex- 
periment using a beam of 130 MeV polarized deuterons 
(equivalent to 65 MeV incident nucleon energy). The ex- 
periment has been performed at KVI in Groningen, em- 
ploying a detector setup covering a large fraction of the 
full breakup phase space. High precision cross sections 
together with vector and tensor analyzing powers have 
been measured in kinematically complete configurations 
by registering energies and angles of the two outgoing 
coincident protons. We have already reported |2| a com- 
parison of the first set of the breakup cross sections with 
realistic TV TV forces and 3NF model predictions, finding 
unambiguously significant effects of the three-nucleon in- 
teraction. In this paper we present an extended breakup 
cross-section data set for 72 kinematical configurations 
(corresponding to a total of about 1200 data points). 
Since we have introduced certain improvements in the 
data analysis procedure, this set partially overlaps with 



the previous data. A second reason for such an over- 
lap is to provide now a fully systematic coverage of the 
phase space, presenting the data on a grid of kinemati- 
cal variables (two proton polar angles, their relative az- 
imuthal angle and the arclength variable). We compare 
our experimental results to theoretical predictions based 
on various approaches. First, we use realistic TV TV poten- 
tials combined with phenomenological 37V interactions. 
Then, we base the predictions on a coupled-channel po- 
tential with the explicit single A-isobar degrees of free- 
dom. Finally, we use also the results of the calculations 
within the ChPT framework at NNLO, with complete 2 TV 
and 3TV dynamics. The comparison is supplemented by 
first global searches of possible regularities in differences 
between the data and theory, determined by inspecting 
the cross sections summed over certain kinematical vari- 
ables. 

There are a few issues which need to be discussed in or- 
der to clarify the details of interpreting the experimental 
results. First of all, we already discussed that the 3TV in- 
teraction is still not completely understood. Recent stud- 
ies of strongly non-local interactions (see 0,E3) aim at 
total removal of the 3NF's. Indeed, the non-locality is 
closely related to the 3NF's 01 and, in principle, this 
can result in an ambiguous separation of 3NF effects and 
off-shell effects. Here we will not discuss this issue fur- 
ther. We only note that our predictions are based on 
several TV TV interactions, some of which are local, some 
are non-local. Nevertheless, they all provide very similar 
predictions, alone and when combined with model 3NF's. 
We also note that the chiral interactions we use are evi- 
dently non-local. It should also be mentioned that all the 
applied formalisms miss two features which are inherently 
present in the experimental data. The first difference is 
the Coulomb interaction: the experiment is performed 
in the deuteron-proton system while all calculations ne- 
glect any long-range forces like the Coulomb interaction. 
It can be argued that the influence of the Coulomb inter- 
action at our energy is (if any) very small. Calculations 
for the elastic scattering cross section at 65 MeV 0, 
indicate an essentially negligible difference for nd and pd 
predictions, even in the cross-section minimum, the most 
sensitive region to study the 3NF effects. The simultane- 
ous treatment of the Coulomb and nuclear forces in the 
Faddeev framework is progressing |5l|. but predictions 
for our breakup data are not available yet. The first in- 
formation suggest, however, that in contrast to the elastic 
scattering case, the Coulomb effects can significantly in- 
fluence the breakup cross sections in certain kinematical 
configurations. Secondly, all the theoretical approaches 
are using a nonrelativistic framework and nonrelativis- 
tic kinematics. Here again we expect the effects induced 
by relativity to be almost negligible. For cross sections 
of the breakup reaction in selected configurations it has 
been shown [52| using relativistic kinematics that the dif- 
ferences between both treatments are minimal at nucleon 
energies below 100 MeV. The remaining problem of ar- 
clength differences does not introduce any noticeable ef- 



fects either - we adopt a projection procedure |4l|, trans- 
forming the theoretical predictions onto the relativistic 
kinematics. Similar conclusions were also reached in [2lJ . 
Ultimately rigorous comparison will be possible only after 
a full relativistic dynamics (boosted potentials) is imple- 
mented for the Faddeev formalism, similarly to the first 
calculations for the 3N bound states 1531 and to the pi- 
oneering Nd elastic scattering study [54J . 

The paper is organized as follows. In Section [n] we 
recall some details of the experiment and of the data 
analysis, emphasizing the refinement introduced since the 
previous report Q. We briefly present in Section ITTll the 
theoretical formalism underlying the calculations based 
on solving the Faddeev equations with the realistic poten- 
tials, with the effective potential obtained in the ChPT 
framework and the coupled-channel approach with the 
explicit A-isobar excitation treatment. Our high preci- 
sion breakup cross-section data are presented and com- 
pared to theoretical predictions in Section llVl We con- 
clude and summarize in Section El 



II. EXPERIMENT AND DATA ANALYSIS 
A. Setup and measurement procedure 

The experiment was performed at the Kernfysisch Ver- 
sneller Instituut (KVI), Groningen, The Netherlands. 
Only the main features of the experimental procedure 
are briefly summarized in the following, the detailed de- 
scription can be found elsewhere [2], l55| . 

The beam of vector and tensor polarized deuterons 
was focused to a spot of approximately 2 mm diameter 
on a liquid hydrogen target of few mm thickness. The 
SALAD (small angle large acceptance detector) de- 
tection system consisted of a three-plane multiwire pro- 
portional chamber (MWPC) and two layers of scintillator 
hodoscope (cf. Fig. QJ. The MWPC was used for pre- 
cise reconstruction of the charged-particle emission an- 
gles. To resolve reconstruction ambiguities for multihit 
events, the MWPC consisted of three active anode planes 
with wires spanned horizontally (x), vertically (y) and 
diagonally (u). The almost point-like reaction region, as 
compared to the target-MWPC distance, allowed for re- 
construction of the polar and azimuthal particle emission 
angles with the overall accuracy of 0.6°. 

The plastic scintillator hodoscope covered the range 
of polar angles between 10° and 35°, and the full range 
of azimuthal angles. It consisted of 24 transmission de- 
tectors (horizontal AE strips) and 24 stopping detectors 
(vertical E slabs), forming together a two-dimensional 
array of 140 elements, with an area of about 60x60 mm 2 
each. The system possessed mirror symmetries with re- 
spect to the horizontal and vertical planes, i.e. it could 
be viewed as composed of four similar sectors, each con- 
sisting of 6 slabs and 6 strips. Strips belonging to one 
sector formed telescopes with slabs of the same sector, 
while they had no overlap with slabs in other sectors. 



4 




FIG. 1: Schematic view of the detection system, presenting 
the relative positions of the MWPC, the two layers of the scin- 
tillation detectors and the scattering chamber containing the 
target. The central beam line transporting the primary beam 
to the target and further downstream to a distant Faraday 
cup is also indicated. For the sake of better view on the E- 
detector wall, 6 AE detectors (one sector) are removed from 
the figure. 



This physical grouping of detectors had a reflection in 
the trigger logic, based on combination of hit multiplic- 
ities within the sectors. Apart from trigger definition, 
information from the telescope array was used for parti- 
cle identification and for determination of their energies. 

The events of interest can be roughly divided into three 
classes. First we distinguish single events, for which only 
one AE-E telescope of the scintillation array has regis- 
tered signals in a proper time window. The other two 
types are coincident events with two particles detected 
in two different telescopes. Among them a distinction 
was made between coincidences of elements belonging to 
the diagonal sectors (candidates for both, elastic scat- 
tering and breakup events) and elements of the adjacent 
sectors or belonging to the same sector (only breakup 
events). These three kinds of triggers were separately 
downscaled, enhancing the coincidence rates, to a level 
acceptable for the data acquisition system. Fine classifi- 
cation of events has been done off-line by incorporating 
the MWPC information - for an extensive description 
see HU. 

For each registered event the information from the 
readout system comprised data from the scintillator ho- 
doscope and from the MWPC. The hodoscope data in- 
cluded times measured with respect to the cyclotron ref- 
erence (rf ) signal and pulse heights for all active detectors 
(strips and slabs). The MWPC information was coded 
into the numbers of the hit wires (more precisely - cen- 
ters and widths of the adjacent groups, i.e. clusters, of 
wires). In addition, several auxiliary pieces of informa- 
tion were stored with each event: the beam polarization 
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state, trigger pattern at various electronic stages, etc. 
Scalers, trigger rates, integrated beam current, pulse gen- 
erator signals for dead-time monitoring, etc., were read 
out every 1 s. 

B. Data analysis 

All basic steps of the data analysis procedure, like 
event selection, energy calibration, determination of de- 
tection efficiencies and cross section normalization, have 
been thoroughly described in Q. The description below 
recalls the main features with an emphasis put on the 
introduced improvements and on additional studies per- 
formed with the aim to reduce experimental uncertainties 
or to control their magnitude with enhanced accuracy. 



1. Selection of events and background subtraction 

The first step of the analysis was an adequate selection 
of the events of interest, i.e. coincident proton - proton 
pairs from the breakup process or, necessary for cross sec- 
tion normalization, deuteron - proton coincidences orig- 
inating from the elastic scattering. To guarantee that 
only the products of the reactions initiated within a sin- 
gle beam burst were selected, a 20-ns-wide time window 
was imposed on the time spectra. Particle identification, 
based on the AE-E technique, proved to be very re- 
liable, providing very good separation between protons 
and dcutcrons in the whole energy range. 

Energy calibration was performed on the basis of data 
collected in special calibration runs with energy de- 
graders of precisely known thicknesses. The positions 
of the peaks corresponding to protons from elastic scat- 
tering which traversed the degraders were compared with 
the results of simulations taking into account all energy 
losses of protons on their paths from the reaction point 
to the detectors. In this way the relation between ADC 
conversion (pulse height) and the energy deposited in the 
E counter was found. The relation between the deposited 
energy and the proton energy at the moment of the reac- 
tion was obtained by analogous simulation of the proton 
energy losses. With all these provisions for each breakup 
event the initial energies of both protons (E\ and £2) 
were determined. 

The coincidence (kinematic) spectra E\ vs. E2 were 
built for each analyzed configuration, defined by polar 
angles 0\, 62 and relative azimuthal angle </>i2 of the two 
emitted protons. The integration limits of A0\ = A62 = 
2° and A(f>i 2 = 10° were used in all experimental inte- 
grations leading to the cross-section results, as well as 
in the studies concerning the performance of the detec- 
tion system. The energies E\ and E 2 of each event were 
transformed into two new variables: D, denoting the dis- 
tance of the (E±,E2) point from the kinematical curve 
in the E\ — E2 plane, and S - the value of the arclength 
along the kinematics. Events in slices along the S axis 




Et (MeV) 

FIG. 2: Ei versus E2 coincidence spectrum of the two pro- 
tons registered at 61 = 20° ± 1°, 2 = 15° ± 1° and <j> 12 = 
100° ± 5°. The solid line shows a 3-body kinematical curve 
calculated for the central values of the experimental angular 
ranges. Inset: determination of the accidental coincidences. 
The spectrum was obtained by choosing one slice along the 
kinematical curve in the E\ vs. Ei spectrum (marked area in 
the main figure) and projecting the events onto the D axis. 
Solid line represents the sum of a linear background function 
(shown with the dashed line) and a Gaussian distribution, fit- 
ted in the range of D between Da and D^, corresponding to 
distances of — 3cr and +3<r from the peak position. 



were projected on the central D axis, as shown in Fig. |21 
In the resulting spectra (inset in Fig. the breakup 
events group themselves in a prominent peak, underlaid 
with only a low background of accidental coincidences. 
As it has been already pointed out in Q, the choice of 
integration limits D a and D^, as well as of the assumed 
background function, is not critical since the contribution 
of accidental coincidences in all analyzed angular config- 
urations was very low (between 2% and 5%). However, to 
treat all configurations in a consistent way, and since all 
the -D-projected distributions have approximately Gaus- 
sian shape, the limits Da and Dj, were chosen at the 
values of — 3cr and +3cr from the maximum of the fitted 
peak. A linear dependence of background between those 
points was assumed. Gaussian shape and linear back- 
ground fitted to a sample distribution are shown in the 
inset of Fig. El 



2. Detection efficiency 

The efficiency in determining the particle-emission an- 
gles, called for simplicity the MWPC efficiency, is a prod- 
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FIG. 3: Global efficiency map for the MWPC, presented as 
a polar plot in the angular coordinates. Only the square-like 
area of the map is meaningful. The range of polar angle 9 
is from 0° to 40° with the binning of 1° and the azimuthal 
angle <f) covers the full range with the binning of 5°. The 
discontinuities of lower efficiency regions are the result of this 
finite binning. 



uct of hardware efficiencies of the MWPC wire planes 
and the efficiency of our procedure of reconstructing 
the angles. Since we accepted only events with ver- 
tical and horizontal wires properly correlated with the 
corresponding E and AE detectors, the ranges of wire 
numbers associated with the individual hodoscope ele- 
ments had to be set wide enough. These correlation ta- 
bles were revised once again (with respect to the proce- 
dure of Ref. inspecting the whole data sample and 
the ranges of wires associated with each hodoscope ele- 
ment have been slightly broadened. In this way, since 
there was practically no uncorrelated noise on the wires, 
no additional background was introduced while the ef- 
ficiency was increased. With these new conditions the 
efficiency has been recalculated in the manner similar to 
the one described in [2] and, in addition, losses due to 
the requirement of correlation between all three planes 
have been determined. In spite of imposing this last re- 
striction with reasonable "safety limits" of ±3 wires in 
u plane, some protons scattered on their way to MWPC 
escaped those limits and were rejected, affecting the to- 
tal efficiency. The final map of the MWPC efficiency is 
presented in Fig. 

Thorough studies of the detection and trigger efficien- 
cies were supplemented with additional tests, performed 
for configurations with the relative azimuthal angle of 
the two protons exceeding 90°. Such configurations can 
be realized by two mutually exclusive classes of events: 
when the two protons were registered in either the adja- 
cent sectors or in two diagonal sectors of the hodoscope. 



FIG. 4: Ratio of event rate of proton-proton coincidences 
registered in two adjacent sectors to the total rate of events 
(adjacent + diagonal) for configuration 6\ — 20° ±1°, 62 = 
15° ± 1°, 0i2 = 140° ± 5°. Points represent experimental 
values obtained by integrating the events within the given 
bin of the arclength S. Solid line shows the level fitted to the 
data, the dashed one the ratio obtained from the simulation. 
Numerical results with statistical errors are given for both, 
experimental and simulated data. 



These two cases corresponded to different trigger signals 
and different downscaling factors, therefore any relative 
inefficiency of the trigger logic and/or of the downscal- 
ing should be reflected by influencing the relative amount 
of events of the two groups. It should be stressed that 
we are sensitive only to the relative efficiency of the two 
trigger classes: since the events originating in the elas- 
tic scattering are always registered in diagonal sectors, 
the trigger/detection efficiency for the diagonal sectors 
cancels out in normalization, cf. Eq. (JIJ. Each of the 
two event groups was analyzed separately and the ratio 
of their rates as a function of the arclength S was con- 
structed. An analogous ratio was calculated for simu- 
lated events. For all the configurations the experimental 
ratio is constant along S and agrees with the result of 
simulations within statistical accuracy of 0.8% or less, 
depending on the configuration (see example in Fig. . 
The above result confirms not only the correct function- 
ing of the trigger but also proper handling of the detec- 
tion efficiencies in the analysis. 

In configurations with the relative azimuthal angle of 
two protons </>i2 not exceeding 90° additional losses of 
acceptance have to be taken into account, due to cases 
when both protons were registered in the same AE or E 
detector (impossible proper particle identification and/or 
energy determination), or in two adjacent E detectors. In 
the latter case, if at least one of the two protons was reg- 
istered close to the edge of the two E detectors, the event 
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FIG. 5: Ratio of the rate of single-sector events to the total 
rate of events (one-sector plus two-sector) for configuration 
0i = 25° ± 1°, 2 = 15° ± 1°, 0i2 = 60° ± 5°. Points show 
the experimental values obtained within 4 MeV S'-bins. The 
difference between the purely geometrical factor (dashed line) 
and the total correction factor (solid line) including the ac- 
ceptance losses for single-sector events, is clearly visible. Nu- 
merical values of the factors are given with their statistical 
uncertainties. 



cannot be distinguished from a so-called cross-over and 
had to be rejected (cf. Sec. lII B 3)l . Obviously, only single- 
sector events are affected by these effects. Therefore, the 
total correction factor is obtained as a product of the 
ratio of events with both protons emitted into a single 
sector to all events collected in the specific configuration, 
and the actual correction, describing the losses within the 
single sector. Experimentally, the total correction can be 
determined by the ratio of the breakup events registered 
in single sectors to all events of the considered configu- 
ration. By means of simulation one can investigate both 
contributions separately. First, in order to find for each 
configuration with ^i 2 < 100° purely geometrical factors 
of probability for single-sector events, an ideal case was 
assumed with no losses due to the detector granularity. 
Then, the events were artificially digitized and analyzed 
in the same way as the experimental ones. It has been 
found that the acceptance losses reduce the geometrical 
factors by up to 13%, depending on the configuration, 
but they are constant as a function of S for the selected 
geometry. In Fig. [S] an example comparing the pure ge- 
ometrical and the total correction factors is shown for 
one configuration. In general the losses increase with de- 
creasing relative azimuthal angle and with the increasing 
difference of the polar angles. It should be noted that 
the errors introduced by applying the correction factors 
are much smaller than the factors themselves. 



3. Cross-over correction 

The procedure of energy calibration faces a problem in 
determining the energy of particles which penetrate from 
one stopping detector to the adjacent one, in the so called 
cross-over events. Simple summing up of the two de- 
posited energies is not completely adequate due to energy 
losses in the foil covering all the detector walls. Addition- 
ally, in a particular situation when the energy deposited 
in one of the E slabs is below the detection threshold, the 
energy information is significantly distorted and, more- 
over, there is even no obvious signature of cross-over. 
Such events are shifted away from the kinematic curve 
and contribute to the background attributed to acciden- 
tal coincidences, which is then subtracted in the way de- 
scribed in Section III B II 

Therefore, a new approach was used, in which all cross- 
over candidates were rejected from the analysis (in a way 
explained below) and their amount was determined with 
the use of Monte-Carlo simulation based on the GEANT4 
package. Narrow regions corresponding to the detec- 
tor borders were defined with the help of high-resolution 
MWPC position coordinates and the particles which en- 
tered those regions and induced signals in two adjacent 
E detectors were discarded. Treating the simulated data 
in the same way it was possible to find the ratio of re- 
jected to registered event numbers for every configura- 
tion, which was used to correct the experimental rates. 
The simulations were performed for elastic scattering and 
for all studied configurations of the breakup reaction. For 
elastic scattering, due to rather high proton energy and 
constrained kinematics the effects are quite large - on 
average about 7% of all the events is biased with the 
cross-over possibility. The corresponding correction fac- 
tors vary strongly, from 4% to 11%, depending on the 
proton polar angle 6 p . Their impact is demonstrated by 
the fact that application of this correction leads to the 
experimental cross-section distribution for elastic scat- 
tering (Fig. with a smoother dependence on 6 p , fol- 
lowing more closely the reference pd data. It is also re- 
flected by a decrease of the x 2 value calculated between 
the two distributions by a factor of about 2 with respect 
to the result obtained for the uncorrected data. Contri- 
butions of the cross-over events, calculated individually 
for each analyzed configuration of the breakup reaction, 
vary between 2% and 5%. The individual cross-over cor- 
rection factors were applied when evaluating the differ- 
ential breakup cross sections, resulting in a decrease of 
their systematic uncertainty. 



4- Cross section normalization 

The breakup cross sections are normalized to the elas- 
tic scattering one, using the measured in parallel rates of 
the elastic scattering events and the available pd elastic 
scattering cross-section data. The differential breakup 
cross section for a chosen angular configuration is thus 
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FIG. 6: Experimental angular distribution of the elastic scat- 
tering events. The absolute normalization of our data (full 
circles) is adjusted to best fit the reference data. Statistical 
errors are smaller than the point size. Empty squares repre- 
sent the reference cross-section values |57f . Solid line shows 
the results of the theoretical calculations obtained with the 
CD Bonn potential and the TM99 3NF model. 



expressed in terms of the elastic scattering cross section 
and both measured coincidence rates: 



d 5 a 



dQid^dS 



(s,ni,n 2 ) 



da, 



dflf 



N br (s,n u n 2 ) 

MhMhAS ' e(fti)e(fi 2 ) 



Atof 



, (1) 



where Nb r is the number of breakup coincidences regis- 
tered at the angles f2i, 2 and projected onto a AS- 
wide arclength bin. Subscripts 1 and 2 refer to the 
first and the second proton registered in coincidence or 
to the proton and the deuteron in the case of elastic 
scattering, f2j = (0j,^>j), with i = 1,2, are the polar 
and azimuthal angles, respectively, and Af2j is the solid 
angle (Af2j = A9iA(f>i sinftj). Products e(f2i)e(02) (or 
e (fif )e (fif )) contain all relevant efficiencies and cor- 



rection factors (cf. Sec. Ill B 311 . N e i is the final number of 
elastic scattering pd coincidences registered at the proton 
angle fif. The elastic scattering cross section ^fj(£lf) 

is taken from j^. The bin width AS was chosen to be 
4 MeV. 

In such an approach we profit from cancellation of 
all factors related to the luminosity, i.e. the integrated 
beam current, the density and the thickness of the tar- 
get. Moreover, since events from both reactions are pro- 



cessed by common electronic and read-out systems, the 
relevant dead-time corrections cancel out in the ratio of 
the registered events. In that way factors which would 
be difficult to determine individually and would induce 
systematic uncertainties are greatly eliminated. 



5. Experimental uncertainties 

A full discussion of the experimental uncertainties has 
been presented in 0] • The statistical accuracy comprises 
the error of the measured number of the breakup coinci- 
dences, as well as statistical uncertainties of all quantities 
used in the cross-section normalization, i.e. the number of 
the elastic scattering events and all efficiencies included 
in Eq. iJTJ. Taking into account the range of the cross- 
section values for our data points, the magnitude of the 
statistical errors varies between 0.5% and 4.0%. 

Non-negligible systematic effects can originate from 
the cross section normalization, from uncertainty of the 
energy calibration parameters, from incomplete cancella- 
tion of polarization effects in the cross section (1.0%) and 
from the procedure of reconstruction of the proton emis- 
sion angles. By introducing the cross-over corrections we 
were able to suppress variations in the experimentally 
obtained elastic scattering distribution and therefore re- 
duce the total uncertainty of the normalization proce- 
dure down to about 2.0%. The uncertainty of the en- 
ergy calibration can result in changing the length of the 
experimental distribution along the S'-curve by at most 
0.7%. The relative cross-section errors resulting from 
such change vary between 0.7% and 2.5%, for central and 
peripheral regions of the measured ^-ranges. This is the 
only systematic uncertainty which changes along the ar- 
clength S in every configuration; all other contributions 
are rather configuration-specific. 

The uncertainty of the reconstructed value of the angle 
is due mainly to finite target thickness, finite size of the 
beam spot on the target, straggling effects and angular 
resolution related to the discreteness of the position in- 
formation delivered by the MWPC. The well-reproduced 
correlation of the proton and deuteron emission angles 
for elastic dp scattering confirms that there is no (at the 
level below 0.3°) systematic shift of the reconstructed po- 
lar angles. The other effects, resulting in smearing out 
of the angular resolution, have been studied by the dedi- 
cated simulation, based again on the GEANT4 package. In 
order to reproduce conditions of the real measurement, 
realistic distribution of the reaction vertices and theo- 
retical angular dependence of the breakup cross section 
were assumed. Straggling effects in materials were intro- 
duced by means of GEANT4 transport routines and, finally, 
the positions of proton trajectories intersecting the wire 
planes were translated to hits on the wires. The same 
reconstruction algorithm as for the real data was applied 
to the simulated events. In this way we were able to com- 
pare the amount of events in each configuration defined 
by proton emission angles, with the number of events in 
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TABLE I: Summary of the relevant experimental cross- 
section uncertainties. Two sample cross-section data points 
(with values close to the minimal and maximal ones mea- 
sured) are selected for presenting individual contributions to 
the systematic uncertainties: 

(1) 0i = 02 = 15°, 012 = 60°, S = 106 MeV, 
dV/dfiid^dS" = 0.078 mb-sr-^MeV -1 , 

(2) 0i = 2 = 25°, 012 = 160°, S = 134 MeV, 
d 5 cr/dfiidfi 2 dS' = 1.57 mb-sr-^MeV" 1 . 

The last column shows all the overall ranges of the relative 
cross section uncertainties. The "total systematic" error is 
obtained by adding the squares of all the contributions. 
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reconstructed angles. It was found that for about 30% of 
events selected on the basis of the reconstructed angles 
the particles were really emitted at angles lying outside 
the chosen angular range. On the other hand, practically 
the same amount of events emitted into the chosen range 
is reconstructed with the values of angles not belonging 
to the considered configuration, and therefore rejected. 
In this way the number of events in "true" and "recon- 
structed" configurations is very well balanced: differences 
are between 0.2% and 1.0% and do not contribute signif- 
icantly to the cross section errors. 

The complete simulation of the breakup process lead 
to the conclusion that the influence of angular resolu- 
tion on the cross section is in fact smaller than what was 
found from the geometrical estimations |2j. Including 
all improvements of the current analysis, the total sys- 
tematic uncertainty is lower by about 1% as compared 
to the previously quoted values. The experimental un- 
certainties relevant for the cross sections presented here 
are summarized in Table [I] The overall ranges of un- 
certainties (last column) are not to be associated with 
particular magnitudes of the cross sections. An obvious 
exception is apparently the statistical accuracy; however, 
since data were collected with different downscaling fac- 
tors and there are certain acceptance losses, this scaling is 
also not straightforward. Therefore, we have selected two 
cross-section points with values close to the minimal and 
maximal of all presented in Sec. IIVI and we display for 
them the individually calculated contributions to their 
uncertainties. One can observe that uncertainties of the 
larger measured cross sections are usually dominated by 



systematic effects, while for the smaller values of cross 
sections the contributions from systematic and statisti- 
cal errors are comparable. 



III. THEORETICAL FORMALISM 
A. Realistic potentials 

The calculation of the cross-section values using re- 
alistic potentials is performed exactly as outlined in our 
previous study , following our standard method for the 
37V continuum. The general overviews of our formulation 
of the 37V scattering problem and of including 3NF into 
the scheme are given in pj and [H^ , respectively. In the 
following a very brief review is presented. 

We use the modern, realistic TV TV potentials AV18 Is], 
charge dependent (CD) Bonn 0,0, and Nijm I and II p. 
Investigating the full 37V system dynamics, we combine 
them with the 27r-exchange TM 3NF, taking its recent 
form [l9| , consistent with chiral symmetry, which will be 
denoted by TM99. The TM99 3NF model contains one 
parameter, At m, used as cut-off to regularize its high- 
momentum behavior. The value of Atm is adjusted for 
each particular combination of the 7V7V force and the 
TM99 3NF to match the value of the 3 H binding en- 
ergy [53. For the four 27V potentials used in the calcu- 
lations the corresponding values of Atm (in units of the 
pion mass m^) are 4.764, 4.469, 4.690 and 4.704, respec- 
tively. 

When the 37V system dynamics is studied with the 
AV18 7V7V potential, we combine it also with the Ur- 
bana IX 3NF [u£ (UIX). To apply it within our frame- 
work it was necessary to transform its configuration- 
space form to momentum space [60| . 

Having the 7V7V and 37V forces, the scattering prob- 
lem in the 37V system is stated in form of a Faddeev-like 
integral equation for an amplitude T: 

T = tP(f> + (l + tG )v£l F (l+P)<l> + tPG T 
+ {l + tG )V$ F (l + P)G T, (2) 

where the initial channel state <j) is composed of a 
deuteron and a momentum eigenstate of the projectile 
nucleon. The TV TV transition operator is denoted by t, the 
free 3TV propagator by Go and P is the sum of a cyclical 
and an anti-cyclical permutation of the three particles. 
The 37V potential V^nf can always be decomposed into 
a sum of three parts: 

V 3NF = V^ F + V$ F + V^ F , (3) 

where the part. V-wf Angles out nucleoli i, on which the 
pion is rescattered. The parts are symmetric under the 
exchange of the two nucleons j and k, with j ^ i 
k. One can see that in Eq. J2J on ly one part, V^ F 
appears explicitly; the others enter via the permutations 
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contained in P. The physical breakup amplitude Uq is 
obtained from T by 



U = (1 + P)T. 



(4) 



Iterating the Faddeev-like equation (J2J and inserting 
the resulting T into Eq. (0} yields the multiple scatter- 
ing series, in which each term contains some number of 
interactions among nucleons via 2N- and 3-ZV-forces with 
free propagation in between. The reaction mechanism is 
thus transparently mirrored. 

We solve Eq. @ using a momentum space partial-wave 
basis pj. To guarantee converged solutions for our case 
of 130 MeV incoming deuterons we take into account all 
partial waves with j max < 6 in the 2N subsystem. This 
gives rise to the maximal number of 142 partial wave 
states in the 3N system for each total 3 AT angular mo- 
mentum J. The convergence has been checked by in- 
specting the results obtained for j max — 6 calculations 
without a 3NF (total number of channels increased to 
194). Finally, the breakup amplitudes Uo have been cal- 
culated for all total angular momenta of the 3A" system 
up to J = 25/2 for any NN interaction, while the in- 
clusion of 3NF's has been carried out for all states up 
to J = 13/2. From this amplitude the cross section is 
obtained in a standard manner [1| . 



B. Chiral Perturbation Theory 

The chiral 2N potential at the next-to-next-to-leading 
order (NNLO) used in the present study is derived from 
the most general effective chiral Lagrangian, based on 
the method of unitary transformation |6l| and using the 
spectral function regularization (SFR) |l3j| . More details 
about the employed regularization schemes and the cor- 
responding cut-offs can be found in 0. Completing 
the 3N system dynamics at NNLO with the naturally 
and consistently arising 3NF contributions is presented 
in 0[. We recall below a few key features. 

The 2N force is obtained by summing up contribu- 
tions from graphs of increasing complexity, accounting 
for, roughly speaking, two kinds of processes: long range 
pion(s) exchanges, where chiral symmetry plays a crucial 
role and short range phenomena, which are effectively 
treated by means of NN contact interactions. The cor- 
responding low energy constants (LEC's) are determined 
from the NN data. The potential is expressed in terms of 
the expansion in powers v of Q/A x , where Q is the soft 
scale, corresponding to the nucleon external momenta 
and the pion mass and A x is the hard scale (around 1 
GeV) associated with the chiral symmetry breaking scale 
or the ultraviolet cut-off (s). For each diagram contribut- 
ing to the potential, the power v can be calculated ac- 
cording to the power counting scheme , see also [(5lJ . 

Up to NNLO, the chiral AW potential can be written 

as 



At the leading order (LO, v = 0) the NN potential V (0) 
is given by the one-pion exchange part (1PE) and two 
contact interactions: 



y(0) = yW y 
' V l7T ' V C 



r<fi) 



(6) 



The leading 1PE term is expressed in terms of standard 
constants: the pion decay constant F n , the pion mass 
m w and the axial- vector nucleon coupling qa ■ In the con- 
tact term two LEC's are introduced, Cs and Ct- The 
next-to-leading order (NLO, v = 2) corrections are due 
to two-pion exchanges (2PE), seven new contact interac- 
tions and a correction to 1PE: 



vl 2) 

v 17T 



(7) 



The leading 2PE term introduces no new parameters (ex- 
cept the SFR cut-off, see below), the contact terms are 
characterized by seven constants Ci,...,CV and in the 
1PE correction term the constant dig can be incorpo- 
rated by renormalization of gA - see 0, |§2 f° r more 
details. Finally, the NNLO (v = 3) corrections are given 
by the subleading 2PE potential and corrections to the 
1PE force: 



v} 3) ■ 



(8) 



there are no new contact terms. The 2PE term contains 
three new LEC's, c\, c 3 and c 4 . The LEC's C's, Ct 
and Ci, . . . , CV, appearing at LO and NLO, are obtained 
by fitting the V NN predictions to the NN data (more 
precisely, to the lowest phase shifts). The three LEC's 
ci, C3 and C4 entering the 2PE contribution at NNLO 
can be determined from the irN scattering data. It has 
been shown 0] that the adopted values lead to a proper 
reproduction of the deuteron properties and of the phase- 
shift analysis results. 

A non-vanishing chiral 3NF arises only at NNLO (in 
the energy-independent formulation) and can be written 
as 
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3NF 



3NF 
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(9) 



yNN = y(0) + y(2) + y(3) 



(5) 



The three terms account for three different topolo- 
gies 0| . The V2^ F describes a simultaneous exchange 
of two pions and incorporates the same LEC's ci 3 4 as 
in the subleading 2PE NN potential V$ . The V^ F 
contribution arises due to a single pion being exchanged 
between a nucleon and a 2A^ contact interaction. The 
contact term contains one parameter, which is usually 
called cd- Finally, the describes a contact inter- 

action of three nucleons and introduces another LEC, 
labeled ce- These two last LEC's are fixed by the re- 
quirement to reproduce the 3 H binding energy and the 
nd doublet scattering length - see [l4| for a detailed dis- 
cussion. 

In representing the chiral potential we use the spec- 
tral function decomposition [16| and we reject the large- 
mass (momentum) fraction of the 2PE via a step Heav- 
iside function with the cut-off parameter A. In the 
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2PE contributions at NLO and NNLO the loop func- 
tions are thus regularized and the corresponding short- 
distance phenomena are shifted into the contact terms 
(the LEC's are appropriately adjusted). This procedure, 
called SFR, possesses numerous advantages over formerly 
implemented dimensional regularization 0, ^| . Its im- 
plementation allowed us also to use LEC's 01,3,4 con- 
sistent with the ttN data, in contrast to the former 
study 14], where the 3N dynamics was described in a 
so-called NNLO* approach, with artificially small values 
of these constants. 

Using the resulting potential, the i-matrix is ob- 
tained via numerical, non-perturbative solution of the 
partial-wave projected Lippmann-Schwinger (LS) equa- 
tion. Since the effective 2N forces are meaningless for 
large momenta, we still have to reject contributions of 
the high-momentum states. In this way we also avoid an 
ultraviolet divergence of the LS equation. The standard 
procedure to accomplish those requirements is to regu- 
larize the potential by multiplying it with a regulator 
function, containin g an additional cut-off parameter A. 
As in other studies [13|, Ha > we use a Gaussian regulator 
function. The T operator of the 3N scattering problem is 
obtained in an identical way as for realistic potentials (cf. 
Sec. Ill Aj l. by solving the Faddeev-like equation (Eq. [2J 
with the chiral 3NF from Eq. (0 . To keep the treatment 
of 2N and 3N interactions consistent, we use an appro- 
priate regulator function with the same cut-off parameter 
A as for the NN potential in the LS equation also for reg- 
ularization of the chiral 3NF. In further calculations the 
observables are generated on the basis of the obtained 
breakup amplitude Uq. 

Our method provides the possibility to estimate uncer- 
tainties of the calculated predictions. We perform calcu- 
lations with a few combinations of the two cut-off pa- 
rameters, [A, A]. The range of predictions obtained for 
reasonable choice of the variation intervals of both cut- 
offs gives an estimate of the theoretical uncertainty. For 
details on how one selects the proper ranges of regular- 
ization cut-off values we refer to [2 ^| and references 
therein. In the present study we use the following pairs 
of the cut-off parameters (values in MeV): 

[A, A] = [450,500], [600,500], [550,600], 

[450,700], [600,700] . (10) 

C. Coupled-channel potential 

A new realistic two-baryon potential coupling NN and 
NA states has been presented in detail in (23 , with sev- 
eral examples of its application to calculate observables 
for the 3N system. The main features of this approach 
are briefly recalled here. 

The dynamics of the 3N system is described with the 
explicit treatment of the A-isobar excitation, considered 
in the relevant energy range as a stable particle. The 
three nucleon channels are coupled to those in which one 



nucleon is excited and forms the A-isobar. Creation of a 
virtual excited state yields an effective 3NF, in parallel to 
other aspects of the dynamics induced by the A-isobar. 

The method for extending a model of NN interaction 
to include A-isobar degrees of freedom has been worked 
out in [6^] and recently thoroughly upgraded [2^] , taking 
the purely nucleonic CD Bonn potential [| as a refer- 
ence. Such a coupled-channel potential is based on the 
exchange of n, p, lo, and a mesons, and in addition to 
the purely nucleonic part includes also contributions from 
the transition between the AW and NA states, from the 
exchange NA — AN potential and from the direct in- 
teraction of the NA states. The force employed here, 
referred to as CD Bonn + A, is as realistic as any of 
the NN force models quoted in Sections |I] and IIII Al 
reproducing the data of the 2N system with a x 2 P er 
degree of freedom of 1.02 |23|. It is purely nucleonic 
in the isospin singlet states; the coupled-channel two- 
baryon extension acts in isospin triplet states only, where 
a few constants of the reference NN force are retuned. 
Prominent contributions of the effective 3NF mediated 
by the A-isobar are of the Fujita-Miyazawa type an d 
of the Illinois ring type [23 • The contributions are based 
on all meson exchanges, i.e. ir, p, a and u> exchanges, 
contained in the coupled-channel potential; the A prop- 
agation is retarded. The arising effective three-nucleon 
force is much richer with respect to the A excitation and 
also has shorter ranged components than standard irre- 
ducible two-pion exchange 3NF's. Furthermore, all its 
components are dynamically consistent with each other 
and with the effective 2N force. However, in addition to 
the A-mediated 3NF an irreducible 3NF covering other 
physics mechanisms is not used. 

The solution of the three-baryon scattering problem is 
based on the AGS equation formulation, using a Cheby- 
shev expansion of the two-baryon transition matrix as 
the interpolation technique . The multichannel tran- 
sition matrix U between two-body channels is obtained 
from: 

U = PG^ 1 +PT a G U 1 (11) 

where T a is the two-baryon transition matrix in three- 
baryon space (the subscript a denoting the pair ^-7 of 
interacting baryons, a ^ (3 ^ 7), Go is the free resolvent 
(E — iO — Ho)^ 1 with the total available energy E and 
free Hamiltonian Hq, and P is the permutation operator 
introduced in Eq. J5J. The transition matrix T a results 
from the full form of the two-baryon potential V a , acting 
between baryons (3 and 7: 

T a = V a + V a Go T a . (12) 

The breakup transition matrix Uq is obtained from U 
according to 

U = (1 + P) Go 1 + (l + P)T a G U. (13) 

The first term on the right side of Eq. (|13|l does not con- 
tribute to the on-shell matrix elements of £/q, needed to 
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calculate breakup observables. The approach described 
here is very similar to the one outlined in Sec. 1111 Al If 
we define the amplitude T as: 

T = T a G U, (14) 

then the integral AGS equation (|ll|l after simple algebra 
becomes identical with the Faddeev-like equation in 
which the 3NF potential is set to zero and the NN t- 
operator is identified with T a . Following the remark be- 
low Eq. H13|l , that equation becomes immediately identi- 
cal with Eq. J3J|. 

Matrix elements of the amplitudes U and T, necessary 
to calculate breakup observables, are found in the partial- 
wave basis. The charge dependence of the two-baryon po- 
tential is treated as described in j23|, yielding the total 
isospin | channels in the state. In the purely nucle- 
onic channels all the states with jmax < 6 in the two- 
baryon system have been taken into account, while for 
the NA channels the applied total angular momentum 
limitation was jmax < 5. In the full three-baryon space 
the states with angular momentum up to J max — 31/2 
were taken into account. For the energy considered here, 
the results are fully convergent with respect to both, j max 
and Jmax limitations, what was tested by checking sev- 
eral predictions obtained with the limits j max — 6 and 

Jmax — 35/2. 

The discussion of the coupled-channel potential ap- 
proach should be closed with a few remarks. The mech- 
anism of explicit A excitation in the three-baryon in- 
teraction has two distinct effects: it yields an effective 
repulsive potential (two-baryon dispersion) and it in- 
duces an effective 3NF. These two contributions usu- 
ally compete [22, , resulting in relatively modest dif- 
ferences when comparing the results of CD Bonn and 
CD Bonn + A predictions. The competition might be 
less pronounced at higher energies. It should be also 
noted that in this method the binding energies of the 3N 
systems are reproduced a bit less perfectly than in the 
other approaches. On the other hand, only in this frame- 
work has a significant development towards including the 
Coulomb interaction into the Faddeev formalism for the 
3A*" continuum been achieved recently [5l| . 

IV. RESULTS 

The main purpose of this paper is a systematic study 
of the quality by which the breakup cross sections can be 
reproduced by theoretical predictions. The investigation 
spans a significant fraction of the breakup reaction phase 
space, the attainable geometries defined by our experi- 
mental conditions. In our methodical approach of scan- 
ning the phase space we present the cross-section data 
for a regular grid of polar and azimuthal angles with a 
constant step in the arclength variable S. Polar angles of 
the two protons B\ and #2 are changed between 15° and 
30° with the step size of 5° and their relative azimuthal 
angle 4>n is analyzed in the range from 40° to 180°, with 



the step size of 20° . We are able to extract data covering 
a denser grid. However, since the changes of the breakup 
cross section are rather smooth, already this coverage 
allows to draw all the important conclusions. For each 
combination of the central values 9\ , O2 and </>i2 the ex- 
perimental data were integrated (cf. Sec. Ill B 1|) within 
the limits of ± 1° for the polar angles and of ± 5° for 
the relative azimuthal angle. The bin size along the kine- 
matic curve S was 4 MeV. Such limits allowed us to reach 
sufficient statistical accuracy while keeping the angle and 
energy integration effects to a minimum, not affecting the 
comparison with the point-geometry theoretical predic- 
tions (see below). 



A. Individual kinematical configurations 

A few geometries on the above defined grid were al- 
ready presented in our previous report |2(. However, 
due to improvements in the analysis procedure discussed 
in Sec. Ill Ul the current results are slightly more pre- 
cise. Therefore, and for the sake of presenting a complete 
picture of data comparison with various theoretical ap- 
proaches, we display in Figs. l7ll4*l cross sections for all 72 
kinematically complete configurations, each figure show- 
ing a collection of 9 geometries (different 0\, 62 pairs) 
for the same value of the relative polar angle $12- The 
data are compared with three sets of theoretical calcu- 
lations, introduced in Sec. 11111 We refer to them by let- 
ters, corresponding to the numbering of the respective 
subsection: the realistic potential approach with model 
3NF's (Se c IhTaI is called "set A", the ChPT predic- 
tions (Sec. IIII Bp are denoted by "set B" and the results 
for the coupled-channel potential with the explicit A- 
isobar treatment (Sec. IIII Cfl are called "set C" . Since 
the predictions of the three sets are often close to each 
other, in order to clearly demonstrate all the details, ev- 
ery figure is composed of two parts. In the upper part 
the data are confronted with the results of calculations 
of set A and set C. The light-shaded bands correspond 
to predictions obtained with only pairwise NN poten- 
tials (AV18, CD Bonn, Nijm I and II), the dark-shaded 
bands show the results when they are combined with the 
27r-exchange TM99 3NF. The dashed lines demonstrate 
the results of calculations with the AV18 potential com- 
bined with the Urbana IX 3NF. The solid lines show the 
predictions obtained with the use of the coupled-channel 
potential CD Bonn + A. In the lower part the same 
data are shown with the predictions obtained at NNLO 
of the ChPT approach. The bands show the ranges of 
the results computed using the different cut-offs, listed 
in Eq. llUfl : the light-shaded bands display the results 
when the calculations were restricted to include only the 
NN force contributions, the dark-shaded bands repre- 
sent the predictions for the full dynamics, with the 3NF 
graphs taken into account. Following the arguments of 
our previous study [2| , we compare the experimental data 
averaged over finite phase-space intervals with the point- 
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FIG. 7: Experimental breakup cross sections in 9 kinematical configurations for the relative azimuthal angle of the two protons 
012 = 40° and for various combinations of the proton polar angles, as indicated in the individual panels. The error bars represent 
statistical uncertainties only. In a few panels the results are scaled with the indicated scaling factors to fit the common vertical 
axis. Upper part: data compared to predictions obtained with the realistic NN potentials only (light-shaded bands), with 
calculations of the pairwise NN forces combined with the TM99 3NF (dark-shaded bands) and of AV18 + Urbana IX (dashed 
lines). The solid line represents the results obtained for calculations within the coupled-channel framework with the CD Bonn 
+ A potential. Lower part: the same data confronted with the predictions obtained within the ChPT approach at NNLO. 
The complete calculations are represented by the dark-shaded bands, while the light-shaded ones demonstrate the results with 
the dynamics constrained to only NN contributions. 
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geometry theoretical predictions calculated at the cen- 
tral values of the ranges of the kinematical variables. It 
has been checked that for all the configurations consid- 
ered here the averaging leads to a slight enhancement of 
the theoretical cross-section values, not exceeding 1.6%, 
equivalent to some extra normalization factor. Since the 
global conclusions are drawn mainly with eliminated in- 
fluence of the data normalization (see further below), this 
simplification does not affect them. 



Figs. 171141 are the basis for the quantitative compar- 
isons of our experimental results with the predictions ob- 
tained in different approaches, as well as between the 
theoretical calculations themselves. 

There is a large number of configurations, concentrated 
mainly (but not exclusively) in the central region of the 
investigated azimuthal angle 0i2 range (Figs. 1101111 top 
panels in Figs. I9I12I14|I . where predictions of all con- 
sidered theoretical approaches are consistent with each 
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FIG. 10: The same as in Fig. [7| but for kinematic configurations with the relative azimuthal angle of two protons <^i2 = 100° 



other over the whole arclength range attainable in our 
experiment. This is particularly true for geometries char- 
acterized by relatively large values of the cross section. 
The bands representing ranges of cross section predicted 
by calculations with different realistic potentials (set A) 
converge practically to a common line, identical with 
the predictions of the coupled-channel potential (set C). 
Similarly, the bands reflecting the computation uncer- 
tainty of the ChPT approach (set B) are also very nar- 



row. Generally, in those configurations the theoretical 
predictions follow very accurately the experimental dis- 
tributions. This confirms the high quality of the pre- 
dictions provided by modern formalisms and simultane- 
ously reflects the precision and accuracy of the experi- 
ment. However, since the predictions with and without 
3NF's are identical, no details of the 37V system dynamics 
can be gained from those data. 

There are, however, regions of the phase space, where 
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the results of the calculations incorporating 3NF con- 
tributions differ substantially from the ones using the 
NN dynamics only. The most pronounced 3NF man- 
ifestations can be observed in the range tpi2 > 120°, 
in the configurations characterized by relatively small 
cross sections (Tigs. lT2*ll4|) . The induced changes concern 
the shape and/or the absolute magnitudes of the cross- 
section distributions. High sensitivity of the predicted 
cross sections to the details of the interaction model ap- 



plied in the calculations makes this region extremely use- 
ful for studying the 3N system dynamics. Regarding the 
results of the realistic potentials approaches (sets A and 
C), one observes that the inclusion of 3NF's usually in- 
creases the predicted cross-section values and that the 
largest effects are introduced by the TM99 3NF model, 
slightly smaller ones for the case of Urbana IX 3NF and 
significantly smaller ones by the explicit A-isobar exci- 
tation. The comparison of the calculated cross sections 
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FIG. 12: The same as in Fig. |7J but for kinematic configurations with the relative azimuthal angle of two protons <^i2 = 140° 



with the data leads to the important conclusion that the 
predictions of the realistic potentials describe the data 
much better when the contributions of the 3NF are taken 
into account. The ChPT results do not reveal that clear 
signal of the 3NF effects. Inclusion of the 3N interaction 
components, if affecting the predictions at all, results in a 
small change of the shape of the cross-section distribution 
along S rather than in modification of its absolute mag- 
nitude. In these geometries the bands representing the 



uncertainties of the ChPT predictions without and with 
the 3NF contributions are relatively wide, they essen- 
tially overlap one another, and generally their agreement 
with the data is satisfactory. One can also observe that 
the calculated cross sections practically coincide with the 
results of the realistic potential approach containing the 
full dynamics, i.e. with the 3NF model included. 

In geometries with small azimuthal angle, 4>i2 < 80° 
(Tigs. 1715)) . one can also identify several cases where the 
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FIG. 13: The same as in Fig. but for kinematic configurations with the relative azimuthal angle of two protons <^i2 = 160° 



contributions of the TM99 or Urbana IX 3NF's modify 
the predictions of the realistic potentials' approach at an 
appreciable level. However, there are few reasons to con- 
sider the situation in this region as qualitatively different 
from that at large <f>i2 values. The coupled-channel cal- 
culations with A-isobar excitation included predict cross- 
section values consistent with those obtained for realistic 
NN potentials without the 3NF contributions. The set 
C predictions even tend to follow the lowest range of the 



set A band (see e.g. configuration Q x = 25°, 8 2 = 25°, l2 
= 40° in Fig.0). Although the ChPT predictions with- 
out and with 3NF contributions still significantly overlap, 
ranges of both kind of predictions are relatively small. 
As for set C, the set B predictions also agree rather well 
with the realistic potentials' results, which do not in- 
clude 37V interaction effects. The comparison of theoret- 
ical results with the data shows noticeable disagreements 
for several geometries. Generally, all approaches, even 
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FIG. 14: The same as in Fig. but for kinematic configurations with the relative azimuthal angle of two protons <^i2 = 180°. 



if their results without and with 3NF contributions are 
almost identical (see e.g. configuration 6\ — 30°, 6 2 = 
15°, 4>i2 = 40° in Fig.[7|), overestimate the data. In ge- 
ometries for which the two kinds of predictions differ, 
this inconsistency is worse for calculations with the 3NF 
contributions (TM99 or Urbana IX) taken into account, 
since adding this piece of dynamics increases the pre- 
dicted cross-section values. 

The discrepancy mentioned above is the largest for the 



configuration Q x = 15°, 8 2 = 15°, 12 = 40° (first panel 
of Fig. 0). All theoretical approaches (although the ef- 
fect is slightly less pronounced for the ChPT predictions) 
deviate from the data by as much as 20%, i.e. far beyond 
the experimental uncertainties. Regarding all configura- 
tions characterized by the smallest analyzed proton polar 
angles, 0\ = 6 2 = 15°, one finds that the disagreement 
between the predictions and the experimental cross sec- 
tion changes systematically, with a strong dependence 
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system. 

Extending the investigation of the ChPT approach be- 
yond the order discussed until now, we have compared 
our data with the predictions including only NN con- 
tributions, obtained at the still higher, N 3 LO order - 
see ^3| for details of the theory involved. Since the 
absence of 3NF contributions makes these calculations 
by virtue incomplete, we include the (incomplete) N 3 LO 
predictions only in global x 2 tests (see below). In Fig. ED 
we present sample comparison of predictions obtained at 
NNLO (dark-shaded) and N 3 LO (hatched) for two ge- 
ometries of the breakup reaction, both based on the NN 
potential only. Comparing calculations based on incom- 
plete dynamics can not be very conclusive, yet we can ob- 
serve the differences between the obtained shapes, what 
may signal the importance of the higher order terms. It is 
also expected that the contributions of the 3NF at N 3 LO 
might be larger than at NNLO. The quantitative compar- 
ison, however, must be postponed until the full dynamics 
of the 3N system is implemented at that order. 



FIG. 15: Comparison of experimental breakup cross sections 
in two kinematic configurations (indicated in the panels) with 
the predictions of ChPT performed at NNLO (dark-shaded 
bands) and at N 3 LO (hatched bands). Both calculations take 
into account only AW force contributions. 



on the relative azimuthal angle: for the small 0x2 val- 
ues the data are overestimated, while for the large 0x2 
they are strongly underestimated. Only for 0x2 around 
100° does the agreement become satisfactory. To a much 
smaller extent this effect is visible also for geometries 
with 9 l = d 2 = 20° and perhaps 6\ = 20°, 6 2 = 15°. It 
should be noted that for all configurations characterized 
by a certain 0±, 9 2 pair the same part of the detector is 
used to extract the data for any 0x2 angle and, more- 
over, the efficiency corrections are tiny in comparison to 
the observed discrepancies, so it is impossible to attribute 
the inconsistencies to any experimental deficiency. 

The presented systematic study, covering a large frac- 
tion of the breakup phase space leads to a rather com- 
plicated picture. Generally, for most of the studied ge- 
ometries the description of data provided by all theo- 
retical approaches is quite satisfactory. There are spe- 
cific regions, where the 3NF effects are pronounced and 
their importance is clearly confirmed by the measured 
cross sections. There are also final state geometries in 
which significant discrepancies between the experiment 
and the theoretical predictions are observed. The pattern 
of disagreement changes as a function of the kinematical 
variables. These findings strongly support the statement 
that only precise measurements in large regions of the 
phase space can provide enough information to judge on 
the quality of the models dealing with the description of 
the breakup observables. Resolving the discrepancies is 
at present not possible; they might be a signal of some 
missing ingredients in the assumed dynamics of the 3iV 



B. Global comparisons 

In order to perform a quantitative comparison of the 
whole bulk of our data with the theoretical predictions 
and to trace possible regularities in (dis-) agreement be- 
tween data and theory, we continued the global tests ini- 
tiated previously , calculating values of x 2 P er degree 
of freedom between the data and individual sets of the- 
oretical predictions. We decided to concentrate on the 
option with a free normalization factor, putting more 
weight onto the shapes of the cross-section distributions 
as a function of S. In this way the conclusions are not 
biased by the absolute normalization uncertainties. In a 
later part of this Section we describe a complementary 
piece of investigation, presenting a comparison of data 
integrated over the arclength variable S with the analo- 
gously treated theoretical predictions. There the exper- 
imental normalization is fully taken into account, while 
the dependence on S (i.e. shape of the distribution) is 
to a large extent neglected. This approach is an exam- 
ple of studying the breakup phase space by inspecting its 
projections onto selected sub-spaces of lower dimensions. 

Focusing our comparison on how the shapes of the 
cross-section distributions are reproduced by different 
theoretical approaches, we have calculated the values of 
X 2 per d.o.f. for all 72 configurations together (a total 
of nearly 1200 cross-section data points) with respect 
to all sets of theoretical predictions. The emphasis on 
the shapes of the distributions is motivated by the ob- 
servation (cf. Sec. IIV A|) that the action of 3NF contri- 
butions is usually equivalent to a small increase of the 
cross-section values in the whole range of S. Therefore, 
if the experimentally determined absolute normalization 
factor would be e.g. slightly too large, the data would 
be artificially shifted towards the predictions including 
full dynamics, leading to erroneous conclusions. This 
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TABLE II: Agreement between the experimental cross sec- 
tions and the theoretical predictions obtained in different ap- 
proaches, quantified in terms of \ 2 P er degree of freedom. 
Major focus in comparing the data with theory is put on the 
shapes of the experimental distributions as explained in the 
text. The quality of the predictions based on only pairwise 
NN interactions (xijv) is compared with X2JV+3JV values ob- 
tained for calculations including genuine 3NF effects, typical 
for the particular approach. Details on the kinds of forces 
used in obtaining \ 2 values at every row are given in the text. 



N N fnrrp 

1 V 1 V Ikjl VjC 


3NP model 


X2N 
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X2N+3N 
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AV18 


TM99 


4.48 


3.80 




Urbana IX 
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CD Bonn 


TM99 


4.04 


3.80 


Nijm I 


TM99 


4.38 


4.38 


Nijm II 


TM99 


4.53 


4.05 
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TM99 
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4.07 






(4.04-4.53) 


(3.80-4.38) 
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ChPT at NNLO 




3.67 


3.96 






(3.36-6.59) 


(3.35-6.64) 


ChPT at N 3 LO 




5.29 






Set C 






Coupled-channel 


A-excitation 


3.83 


3.63 



method also allows us to eliminate the small influence of 
averaging, inherently present in the data and omitted in 
the theoretical predictions (averaging does not affect the 
shapes of the cross-section distributions presented here) . 
To eliminate the influence of the absolute experimental 
normalization, the data were renormalized in each con- 
figuration by a constant factor (limited to the range 0.9 
to 1.1), to best fit the particular theoretical distribution. 
In this way the quality with which a given set of theo- 
retical predictions reproduces all the data is quantified 
by a single number. In particular, for each combination 
of forces we can compare two x 2 values: x\n obtained 
for predictions based on pairwise NN interaction only 
and X2N+3N f° r the calculations including 3NF contri- 
butions. It should be noted that in the x 2 analysis only 
statistical uncertainties were taken into account, there- 
fore values exceeding 1 can be expected. Investigating 
influences of the 3NF effects we concentrate rather on 
the relative change from \\n *° X-2N+3N ano - n °t on the 
absolute x 2 value. The same argument holds when com- 
paring predictions of different forces with respect to the 
quality with which they describe the experimental data. 

The results of the x 2 analysis for all considered theo- 
retical approaches are shown in Table [H] First, the two 
kinds (without and with 3NF contributions included in 
the theory) of x 2 values are shown for 4 realistic NN 
potentials and their combination with the TM99 3NF. 
The second row for the AV18 potential gives X2N+3N f° r 



this force combined with the Urbana IX 3NF. Exclud- 
ing this last combination, we define a "mean realistic" 
prediction as a set of cross-section values given at each 
point (Oi,02,4>i2,S) as a mean between the minimum and 
maximum cross section predicted by the 4 realistic forces 
(or their combination with TM99 3NF) at this kinemat- 
ical point. The x 2 values with respect to this mean re- 
alistic prediction are also shown in Table [n] with the 
ranges, equal to the corresponding extreme x 2 values, 
repeated in the next row. These values are to be com- 
pared with the ones obtained for the ChPT calculations. 
In this case only the x 2 for the "mean" set is quoted. It 
is obtained analogously as in the realistic potentials case, 
as the central value between extremes predicted with 5 
combinations of cut-offs. Also the ranges of x 2 values for 
predictions at NNLO are shown for comparison. In the 
case of N 3 LO calculations we quote only the central x 2 > 
reminding the reader that the bands based on the NN 
forces only tend to be wider at this computational order 
and therefore also the accuracy of the predictions would 
behave accordingly. The last row of Table [Q] presents 
the results of x 2 analysis performed with respect to the 
coupled-channel potential. A small difference in xIat va l~ 
ues for CD Bonn between set A and set C calculations 
is due to the different treatment of the charge depen- 
dence. The X2N+3N denotes here the value obtained for 
the calculations including all A-isobar excitation effects. 

Comparing the numbers presented in Table [H] one ob- 
serves that combining any of the realistic potentials with 
a 3NF model improves the description of our data, de- 
creasing the x 2 by about 10%. The effect varies for dif- 
ferent NN potentials of the set A calculations so that the 
ranges of \\n an d x\n+3N va l ues overlap. Nevertheless, 
a systematic shift of the predictions towards the data is 
visible for calculations including the full dynamics. Pre- 
dictions obtained within the ChPT framework do not 
allow for such a conclusion. Although the quality of the 
data description for the set B calculations is very similar 
to that of the realistic potentials approach, the ranges of 
X 2 values obtained with and without 3NF contributions 
are very wide and overlap completely. The central- value 
predictions reveal even a slight worsening of the data de- 
scription induced by including the 3NF effects. However, 
this observation is rather far fetched in view of large the- 
oretical uncertainties and, moreover, taking the center of 
the relatively wide bands is not relevant for tracing de- 
tails of the shapes of the distributions. The predictions of 
NN force alone at N 3 LO show a much poorer agreement 
with the data than those of NNLO. Since the total 3NF 
contributions up to that order are expected to be larger 
than at NNLO, this effect might be easily compensated 
for by the complete dynamics included in the formal- 
ism. The set C calculations, with the coupled-channel 
potential and explicit A-isobar degrees of freedom, lead 
in general to the smallest values of x 2 - The effects of the 
A excitation are rather modest, but with their inclusion 
the predictions are moved by a few percent closer to the 
data. This conclusion is, however, biased by the absence 
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of any estimate of the theoretical uncertainties. 

As mentioned in Sec. IIV Al the largest disagreements 
between the data and the theoretical predictions are ob- 
served for configurations with the smallest polar angles, 
9i = 02 = 15°. In order to eliminate their possible dom- 
inant impact on the x 2 analysis, we have recalculated 
all the values of Table [D] excluding this piece of data (8 
configuration out of 72). It has been found that the x\n 
and X2N+3N va lues obtained in that way were decreased 
only by about 5%, but the overall picture was preserved 
and thus all the above conclusions are valid also for that 
limited data sample. Below we present less global com- 
parisons based on this subset of the data. 

To search for possible regularities in changes of the 
quality of the data description by the models, a less global 
treatment is obviously needed. Firstly, we studied the 
consistency between the data and theoretical predictions 
in various regions of the phase space, inspecting the de- 
pendence of x 2 on the relative azimuthal angle <f>i2 of the 
two protons. Values of x 2 have been calculated as in the 
global comparison case, but for the groups of configura- 
tions characterized by the same </>i2 value, i.e. separately 
for each group presented in Figs. 171141 The results are 
shown in Fig. 1161 by three different symbols represent- 
ing the three calculation sets and the x\n an d X2N+3N 
shown by open and full symbols, respectively. To sim- 
plify the picture the realistic potentials are represented 
by the CD Bonn force (with and without TM99 3NF) and 
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FIG. 17: Ratios of x 2 values for calculations without and with 
3NF contributions. The cross-section data with the experi- 
mental absolute normalization were used in computing the 
X2N an d xliv+3Af values for groups of kinematical configura- 
tions with the same <j)i2 angle. The results of data comparison 
with calculations of sets A, C and B are shown by dots, tri- 
angles and squares, respectively. The lines are only to guide 
the eye. 



the ChPT approach results are shown for the mean set. 
One observes that for 4>i2 < 90° there is practically no 
effect of including the TM99 3NF into the calculations. 
For larger relative azimuthal angles the description of 
data is significantly improved by employing the full dy- 
namics. The coupled-channel calculations predict much 
smaller effects due to A excitation, what is the result of 
a compensation mechanism, cf. Sec. lIIlCl Moreover, the 
quality of description depends only very weakly on (^12 . 
The ChPT predictions do not show any large effects of 
3NF contributions either. Only for the largest azimuthal 
angles, 4>i2 > 120°, in contrast to the CD Bonn results, 
the full dynamics reproduces the data worse than calcu- 
lations with NN interaction terms only. 

Pursuing the study of the x 2 dependence on the rela- 
tive azimuthal angle, we have also checked the changes 
of quality of data description by calculations with and 
without 3NF contributions for the cross sections with 
absolute experimental normalization applied to the data. 
We present the results as the ratio of \\ N to x\n+zn 
in order to magnify the influence of the 3NF effects. In 
Fig. El the ratios for the same theoretical sets as for the 
free normalization case are shown with the same symbols. 
One finds again that the consistency between the predic- 
tions of the CD Bonn potential and the data is improved 
by adding TM99 3NF in configurations with relatively 
large <pi2 angles (ratio above 1). On the other hand, for 
012 < 90° including the 3NF into the calculations moves 
the results away from the data. Astonishingly, the mag- 
nitude of the relative change is almost the same in both 
directions, described by a factor of about 2. For the 
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ChPT calculations no effect is present, the ratio stays 
close to 1 for all values of 12 . The behavior revealed 
by the set A predictions is qualitatively confirmed by the 
coupled-channel calculations, however the amplitude of 
the changes induced by including the A-excitation con- 
tributions is smaller. 

A great advantage of an experiment with the position 
sensitive detector covering a significant part of the phase 
space is the opportunity to study dependences of the ob- 
servables (here the differential cross section) on all inde- 
pendent kinematic variables. However, inspecting the re- 
sults in many-dimensional space is difficult and the com- 
parisons with the theoretical predictions might miss the 
regularities. One possible solution to reduce the complex- 
ity of the problem and still be able to use all the data is to 
select a small number (e.g. 1 or 2) of variables and to in- 
tegrate the observable over the others. Integration of the 
experimental data is usually quite straightforward - it is 
accomplished by summing up of events which fulfill the 
required conditions. But these experimental conditions 
(acceptances, thresholds, granularity, etc.) impose limi- 
tations, which make the procedure of integration for the 
theoretical predictions very complicated and the compar- 
isons might be jeopardized by introducing uncontrollable 
systematic errors. A possible method to resolve the prob- 
lem of comparing the integrated experimental and theo- 
retical observables has been suggested in j6j|. It allows 
to effectively integrate the calculated observables over all 
but one variables, with all experimental constraints taken 
into account. In the case of cross sections, however, such 
an approach leads to numerical values which are hard to 
interpret physically. Influences of the physical changes 
(due to reaction dynamics) of the observable are merged 
with the acceptance functions and the comparisons are 
meaningful only for "integrated physical values" . There- 
fore, in our first attempts to investigate regularities in the 
breakup phase space we rather employ a simpler method, 
deconvoluting the acceptances from the experimental re- 
sults and comparing the integrated cross sections with 
the accordingly summed theoretical predictions. In this 
way we end up with "objective" cross-section values, with 
direct physical interpretation. 

The cross-section results for the individual configura- 
tions shown in Figs. I7IT41 suggest a possible correlation 
between the polar and azimuthal angles with respect to 
the quality of the agreement between the data and pre- 
dictions. Therefore, we studied the cross-section depen- 
dences on the proton emission angles with the experi- 
mental data integrated only over S. The energy thresh- 
old of the detection system introduces an inherent in- 
fluence of the instrumental acceptance onto the result, 
but it is easy to reproduce without any detailed knowl- 
edge of other features of the apparatus. To guarantee an 
exact equivalence of the low-energy cut-off condition for 
the experimental and the predicted results, a threshold 
of 25 MeV (higher than the hardware level) was applied 
for both proton energies and only the range of S lim- 
ited by this requirement was included in the integration. 
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FIG. 18: Differential, integrated over S, cross-section values, 
presented as functions of the relative azimuthal angle cj>i2, for 
several pairs of the proton polar angles 6\ and #2 (indicated 
in the panels). The integration limits are defined by setting 
the proton energy threshold at 25 MeV. The data points are 
compared with the results of calculations with the CD Bonn 
potential (dashed lines) and with the CD Bonn + TM99 3NF 
combination of forces (solid lines). 



The results for 6 pairs of the proton polar angles as func- 
tions of the relative azimuthal angle are shown in Fig. (T51 
The integrated experimental cross sections are compared 
to the correspondingly integrated theoretical predictions 
based on the CD Bonn potential only (dashed lines) and 
with the TM99 3NF included in the calculations (solid 
lines). It has been checked that the results for other re- 
alistic forces are almost indistinguishable from the ones 
presented in Fig. The tendency, already visible in 
the cross-section plots for individual configurations, can 
be better traced here. For 4>i2 below 90° the theoretical 
predictions overestimate the data and the discrepancy 
rises with decreasing 4>\2 values. In the central region of 
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the analyzed 4>i2 range the agreement between the data 
and the theoretical curves is the best. With further in- 
crease of the 0i2 angle the theoretical predictions start 
to underestimate the data. This discrepancy is, however, 
reduced (in various fractions) by including the 3NF into 
the calculations. On the contrary, for the small 4>i2 an- 
gles the effects of 3NF inclusion increase the discrepancies 
between the predictions and the data. It can be finally 
stated that in all cases studied in Fig. ^] the slope of 
the data, though qualitatively reproduced, is not exactly 
matched by the theoretical cross-section values and that 
overall rather small 3NF contributions do not change the 
global picture. Again, this has to be attributed to some 
still unresolved deficiencies of the models of the 37V sys- 
tem dynamics. 

Recapitulating all the results, from both the individ- 
ual configurations of Figs. 171 FH and the global tests from 
this section, we can conclude, that the present day mod- 
els of the 3iV system dynamics reproduce the majority 
of the data with satisfactory precision. In many cases 
in which the predicted effects due to 3NF's are non- 
negligible, their inclusion tends to improve the agreement 
with the data. However, thanks to the applied experi- 
mental technique of covering a significant fraction of the 
breakup phase space with a highly symmetric detection 
system, it has been shown that there are also systematic 
regularities in the discrepancies between the measured 
cross sections and the predictions of all the theoretical 
approaches. Since the systematic factors of our measure- 
ment are common to all configurations, the established 
trends cannot be attributed to systematic experimental 
uncertainties and therefore hint at missing ingredients 
of the nuclear Hamiltonian models. It should also be 
stressed that additional complete and precise data sets, 
at other energies and in even larger phase space regions, 
are needed to study details of the interactions in the few- 
nucleon system. 



V. SUMMARY AND CONCLUSIONS 

A measurement of the deuteron-proton 1 H(d,pp)n 
breakup cross sections using a 130 MeV deuteron beam 
was performed for a large part of the available phase 
space. In this paper high precision, five-fold differential 
cross-section data for 72 kinematically complete configu- 
rations (total of nearly 1200 cross-section data points) at 
different angular combinations of the two outgoing pro- 
tons are presented. We discuss first examples of a global 
analysis of the data, which is trying to establish possible 
regularities of the (dis-)agreement between the experi- 
mental data and different theoretical approaches. 

We compare the measured cross sections to theoretical 
predictions treating the full dynamics of the 3 N system in 
three different ways: employing the realistic NN poten- 
tials AV18, CD Bonn, Nijm I and Nijm II and including 
the 3NF effects by combining them with the TM99 3NF 
model (for AVI8 also Urbana IX 3NF), obtaining the 



nuclear effective potential in the ChPT approach with 
the calculations performed at NNLO order with NN and 
3N contributions (mentioning also the pure JVJV results 
obtained at N 3 LO), and using the coupled-channel tech- 
nique of explicit inclusion of a single A-isobar degrees 
of freedom, resulting in a modified form of the realistic 
CD Bonn force and its combination with all single A ex- 
citation effects in the three-baryon system. The three 
approaches match equally well the properties of the JVJV 
system. When only the 27V dynamical sector is used, 
their predictions for the breakup cross sections are es- 
sentially equivalent. 

The magnitude of the predicted 3NF effects depends on 
the approach. In the case of the coupled-channel poten- 
tial the influences of the A-isobar excitation are generally 
rather small. This is due to a competition of effects in- 
duced by two mechanisms, the two-baryon dispersion and 
the effective 3NF. ChPT predictions, considered in terms 
of ranges of the cross-section values computed with dif- 
ferent cut-off parameters, reveal also rather weak contri- 
butions of the 3NF effects. They are usually smaller than 
the residual dependence on the cut-offs. Also the calcu- 
lations at N 3 LO (presently available with JVJV contribu- 
tions only) are characterized by a broad range for the 
predicted cross-section values. The deviations of NNLO 
and N 3 LO predictions depend on the final state geome- 
try. The largest sensitivity to the 3NF is found in the 
realistic potentials approach. But even for the realistic 
forces there are several final state geometries in which the 
3NF effects are practically negligible. Generally, in such 
cases the cross-section data are in good agreement with 
the theoretical predictions. However, in many analyzed 
configurations the effects of including the 3NF are not 
negligible. Taking into account the 3NF contributions in 
the calculations leads to an increase of the cross-section 
values. This effect is slightly less pronounced for com- 
bining the Urbana IX 3NF with the AVI8 NN potential 
than for the TM99 force combined with any of the four 
considered NN potentials. 

A global analysis, focused on the shapes of the cross- 
section distributions as functions of the arclength vari- 
able S, shows that the agreement between the experimen- 
tal data and the theoretical predictions improves when 
the 3NF contributions are taken into account. This con- 
clusion is valid for all combinations of realistic NN po- 
tentials with the model 3NF's. While for the ChPT pre- 
dictions no conclusions can be drawn due to essentially 
overlapping ranges of predictions without and with 3NF 
included, the coupled-channel calculations also reveal a 
slight improvement in the description of the data when 
the single A-isobar excitation effects are incorporated. 

There is, however, a number of configurations in which 
the cross-section data are not correctly reproduced by 
any calculation. The effect depends on the relative az- 
imuthal angle <f>i2 of the two protons: for small values 
the data are overestimated by the predictions, the agree- 
ment becomes good in the central range of the analyzed 
012 and at the largest angles the discrepancy is reversed. 
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Thanks to the highly symmetric form of our detection 
system, which allows to reduce and to carefully control 
and test for systematic uncertainties, we basically ex- 
clude the possibility of attributing this inconsistency to 
any experimental deficiency. 

The regularities of disagreements have been studied in 
more detail using a global analysis, in which we have 
concentrated on both, the shapes of the distributions as 
well as on the absolute values of the cross sections. It 
has been established that for configurations with large 
values of the <f>i 2 angle the agreement between the exper- 
imental data and the theoretical predictions is improved 
when the 3NF contributions are taken into account. On 
the contrary, for <j>i 2 < 100° the 3NF effects move the 
predictions away from the data. This conclusion is valid 
for all combinations of realistic NN potentials with the 
model 3NF's. It also holds for the predictions obtained 
in the coupled-channel approach, but with reduced size 
of the effects induced by the A-excitation. The ChPT 
calculations predict essentially no sensitivity to the 3NF 
influences along (pi 2 . 

We have confirmed that sizable influences of 3NF's 
are visible in the breakup cross sections at the energy 
of our measurement. Comparison of the agreement be- 
tween the experimental data and the predicted cross- 
section values is presented by first examples of analyz- 
ing a multi-dimensional breakup data set by inspecting 
its projections onto selected sub-spaces of lower dimen- 
sions. Since the advantages of the experimental method 
reduce strongly the impact of systematic errors, the es- 
tablished discrepancies might be considered as signal of 
some missing ingredients in the assumed dynamics of the 
3iV system. Determining regularities in the disagree- 
ments might eventually help to identify shortcomings of 
the nuclear force models. However, it cannot be ruled 
out that the discrepancies result from Coulomb effects, 
which arc ignored in all theoretical approaches presently 
used. In view of the recent advances in including long 
range forces into the Faddeev formalism this persisting 
question might soon be quantitatively addressed. 

Our present study clearly shows the importance of 
complete, precise experiments, determining various ob- 



servables of the breakup reaction. Implementing sym- 
metric detection systems covering large fractions of the 
phase space allows to attain rich data sets, for which sys- 
tematic uncertainties are strongly suppressed and well 
controlled. Results of such experiments are the basis 
for systematic comparisons with different theoretical ap- 
proaches. They provide stringent criteria for verification 
of the models of the nuclear Hamiltonian, assumed in 
formulating the three-nucleon scattering problem. Fur- 
ther improvements of the theoretical models, which are 
also a basis for exact calculations in many-body systems, 
require still more experimental input. We are going to 
supplement the cross-section results with polarization ob- 
servables, with which we will be able to provide more de- 
tailed information, and hopefully pin down the discrep- 
ancies found here. More data sets are needed, acquired 
at energies lower and higher than in our study, with pro- 
ton and deuteron beams, as well as covering still larger 
fractions of the phase space. 
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